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I. INTRODUCTION 

The precise shape of the quark-gluon vertex plays an 
essential role in the understanding of many important 
properties of the infrared behavior of QCD [ll-y|- This 
vertex is in particular fundamental in the analysis of the 
origin of dynamical chiral symmetry breaking in the con¬ 
text of truncated Dyson-Schwinger (DS) equations. In¬ 
deed, the DS equation for the quark self-energy requires 
the knowledge of this vertex, together with the gluon two- 
point vertex function. Even if this kind of object depends 
on the gauge-fixing, one expects that the gauge-invariant 
physical content of the theory is encoded on its behavior. 
For example, the meson spectrum can be studied by in¬ 
jecting such vertices and two-point functions in appropri¬ 
ate Bethe-Salpeter equations, see for example For 

these reasons, the determination of the quark-gluon ver¬ 
tex has been a topic of intensive work in the past. It has 
been first studied in the perturbative domain. A general 
parametrization based on its symmetries and a one-loop 
calculation of the abelian part in Feynman gauge was per¬ 
formed in 0. Some years later, the complete one loop 
calculation in an arbitrary linear covariant gauge and in 
arbitrary dimension was performed in [l^ (where pre¬ 
vious partial results are reviewed also). More recently 
the complete 2-loop off-shell calculation has been per¬ 
formed in d = 4 These perturbative studies lead to 
a satisfactory description of this vertex at high momenta. 
However, the infrared sector, which is of direct interest 
for the understanding of chiral symmetry breaking, is not 
accessible to perturbation theory. To cope with problem, 
mainly two methods have been used: lattice simulations 
and DS equations. 

In both of these approaches, the most convenient gauge 
is the Landau gauge. From the lattice simulation side. 


it is well-adapted because the gauge-fixing operation can 
be expressed in terms of an extremization, which can be 
implemented efficiently on the lattice. From the analytic 
side, this gauge has several advantages. First, it is a 
linear gauge condition that manifestly preserves Lorentz 
invariance. Second, among all linear covariant gauges, 
it is the only one that includes the symmetries of the 
Curci-Ferrari gauge [T3 - [l^ . This implies, in particular, 
the existence of some non-renormalization theorems, that 
simplify considerably the algebraic work [TtI - I^ . Third, 
in this gauge, the gluon propagator is transverse. Conse¬ 
quently, only the vertices where the gluon lines are mul¬ 
tiplied by transverse projectors contribute to the correla¬ 
tion functions. Accordingly, all consideration below will 
be done in this gauge fixing. 

There have been many analysis of the infrared behav¬ 
ior of the quark-gluon vertex based in the study of trun¬ 
cated DS equations (for a review, see (1| ; for more recent 
work on the subject, see (2ll - [^ ). In parallel the inter¬ 
mediate and infrared regime of the vertex has been stud¬ 
ied by lattice simulations. Many tensorial components 
have been analyzed for some particular configura¬ 

tion of momenta. Moreover, for the tensorial component 
that is already present at bare level, which is expected to 
be “dominant”, more gen eral configurations of momenta 
have been considered (^ . 

In continuum analytical approaches, the gauge fixing is 
generally implemented through the Faddeev-Popov con¬ 
struction, which is a perfectly valid procedure in the 
ultraviolet regime, where standard perturbation theory 
works well and where the Faddeev-Popov construction is 
well justified. However, for momenta smaller or of the 
order of 1 GeV (which is the typical energy scale of the 
strong interaction) this procedure is not fully justified. 
The standard lore is that it becomes inappropriate be- 
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cause the coupling grows and perturbation theory even¬ 
tually becomes a very bad approximation scheme. This 
argument must be tempered because the real expansion 
parameter of QCD is of the order of A^as/(47r) and be¬ 
cause lattice simulations, at least in some renormaliza¬ 
tion schemes, indicate that this parameter remains mod¬ 
erate (see, for example [13 )■ Accordingly it is reasonable 
to believe that some sort of perturbation theory should 
give results that are in qualitative agreement with exper¬ 
iments or lattice simulations in that regime. There is, 
however, another reason why standard perturbation the¬ 
ory should become inappropriate in the infrared regime: 
the Faddeev-Popov procedure itself becomes question¬ 
able. This originates in the existence of Gribov copies 
[ 2 ^ in covariant gauge-fixings of non-abelian gauge theo¬ 
ries that are ignored in the Faddeev-Popov construction. 
These copies play no significant role in the ultraviolet 
regime but a systematic treatment of the associated am¬ 
biguity remains an open problem in the infrared. 

The most developed procedure which aims at taking 
into account the effects of Gribov copies is probably the 
Gribov-Zwanziger (GZ) construction [29I - I33 - The pur¬ 
pose of this approach is to restrict the functional integral 
to the first Gribov region where the Faddeev-Popov op¬ 
erator is positive definite. Gribov showed that in this 
region there are no copies that are an infinitesimal trans¬ 
formation of one another. Unfortunately it was shown 
later that there are finite gau ge t ransformations that pre¬ 
serve the Landau condition |33|. Gonsequently the GZ 
procedure does not remove all copies. Moreover, the im¬ 
plementation of the restriction of the functional integral 
to the first Gribov region is not completely justified from 
first principles. 

Despite these difficulties, first studies of the GZ pre¬ 
dicted that the gluon propagator would vanish at zero 
momentum, instead of diverging as in the bare theory. 
A similar behavior was found also as solutions of trun¬ 
cated DS equations P, IMlil and in Non-perturbative 
Renormalization Group (NPRG) calculations (for 

a general review on this formalism, see M)- This kind 
of solution is called “scaling solution”. The suppression 
of the gluon propagator was rapidly confirmed by lattice 
simulations. However, when lattice simulations became 
more precise they showed that the gluon propagator in 
Landau gauge did not behave precisely as the “scaling 
solution” that DS had suggested, neither in quenched 
lattice simulations [dol - li^ nor in the presence of dynam¬ 
ical quarks [13, In fact, lattice simulations favored 
a gluon propagator which tends to a finite value at low 
momentum in d = 4. This massive behavior (sometimes 
called “decoupling” solution) gives a propagator which 
is suppressed with respect to the bare one but not as 
strongly as in the “scaling solution” of DS equations. 
In parallel such decoupling solutions were also obtained 
from DS and NPRG equations [s^, l49l - [^ , as well as in 
a “refined” version of the GZ construction . 

Nowadays, all approaches finally converge toward 
a massive-like gluon propagator in Yang-Mills theory. 


Moreover, in what concerns the ghost sector, there is 
also a consensus in favor of a massless behavior, with a 
slight enhancement at low energy but with the tree level 
power-law. Finally, we note that DS and NPRG equa¬ 
tions favor a coupling constant bounded in the infrared 
at not too large values, see for instance (55| . 

All these observations led two of us to test wether a 
simple perturbative analysis where a massive gluon prop¬ 
agator was introduced at the bare level could repro¬ 
duce quantitatively the lattice data for several correlation 
functions. This lagrangian is invariant under a modified 
Becchi-Rouet-Stora-Tuytin (BRST) symmetry that en¬ 
sures the renormalizability of the model (Sa - I^ . How¬ 
ever, because of the gluon mass, the associated trans¬ 
formation is not nilpotent. As a consequence, the text¬ 
book construction of a physical space based on the co¬ 
homology of the BRST charge does not apply. It must be 
stressed that this difficulty is common to all methods that 
go beyond the standard perturbative FP method. In par¬ 
ticular, the standard definition of the physical space does 
not apply to the GZ model but even stronger, does not 
apply to the procedure implemented on Landau-gauge 
lattice simulations. In fact the standard perturbative def¬ 
inition based on the cohomology of the BRST operator 
includes in the physical space the transverse gluons. But 
this is not satisfactory for two reasons. First, on phys¬ 
ical grounds, we know that gluons are confined and are 
not part of the physical spectrum. Second, if trasverse 
gluons were in the physical space, the spectral density as¬ 
sociated with their two-point function should be positive 
and lattice simulations deary show positivity violations 

[liliS- 

All these considerations show that a proper definition 
of a physical space is an open and important problem. 
In this sense, the lack of a nilpotent BRST charge in 
the Curci-Ferrari should not prevent us from studying 
it in relation with the lattice simulations of correlation 
functions. It is an open question to know if one can define 
a physical space for the Curcci-Ferrari model where all 
states would have a positive norm. 

In the present article we exploit this massive extension 
of the Faddeev-Popov Lagrangian in order to study the 
quark-gluon vertex in Landau gauge. Our aim is twofold. 
First, we want to compare our findings with lattice simu¬ 
lation in order to have a further test of the ability of the 
Curci-Ferrari model to reproduce lattice simulations. As 
we show below, we are able to compute the full vertex 
at one loop in any dimension, for all tensorial structures 
and for any configuration of momenta. The comparison 
of this one loop result is in reasonable agreement with 
all available lattice data (in many cases the agreement 
is very good). In particular, we discuss in detail the re¬ 
sults of a certain scalar function (A 2 in the terminology of 
M) where a significant discrepancy has been observed 
between lattice and DS results [13, We discuss this 
discrepancy in detail giving a very simple explanation of 
its origin. Second (and more importantly) as explained 
before, this vertex is very important in various physi- 
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cal applications of DS and Bethe-Salpeter equations. We 
expect that the present result will be useful in future 
physical applications in these equations. 

We want to stress here that, as it stands, the sim¬ 
ple one-loop calculation presented in this article does 
not encode a spontaneous breaking of the chiral sym¬ 
metry. However, the renormalization-group flow leads to 
a strong enhancement of the quark mass in the infrared, 
which mimicks, for all practical purposes, the behavior 
induced by this symmetry breaking. This is a key ingre¬ 
dient for reproducing the infrared behavior of the quark- 
gluon vertex. 

The rest of the article is organized as follows. In Sec.ITTl 
we describe in more details our model and shortly re¬ 
view the systematic comparison of one-loop computa¬ 
tions with lattice data. We give some details on the 
one-loop calculations in Sec. IIIII and finally present our 
results in Sec. El 


II. THE MODEL 

In this article, we pursue our systematic analysis of the 
correlation functions of QCD in the Landau gauge, in a 
scheme where perturbative calculations can be performed 
in a controlled way for all energies. As explained in the 
introduction we consider the (euclidean) Curci-Ferrari la- 
grangian 0 : 

1 77 ? ^ 

Za =-F%K^ + 

( 1 ) 

where is the gauge field, c“ and c“ are ghost and 
antighost fields, /i“ is a Lagrange multiplier that enforces 
the Landau gauge condition = 0, g is the coupling 

constant and 

FI, = d^Al - d.Al + gr^^AlAt, 

[D^cT = d^c<^ Fgr’^'^Ay. 

It is important to recall that mass is introduced at the 
level of the gauge-fixed action. As a consequence, the 
bare propagators decrease at large momenta as 1 /p^ and 
the mass induces a soft breaking of the BRST symme¬ 
try. As a consequence, all the ultraviolet properties of 
the standard Faddeev-Popov action are preserved to all 
orders in perturbation theory, in particular renormaliz- 
ability and the behavior of correlation functions. 

One of the main consequences of the mass term is that 
the low energy properties of the system are more regular 
than within the usual Faddeev-Popov approach. In par¬ 
ticular, it is possible to find renormalization schemes in 
which the coupling constant remains finite down to the 
deep infrared [^, at odds with the findings of standard 
perturbation theory which presents a divergence of the 
coupling constant at some energy of the order of 1 GeV. 
The absence of such a Landau pole opens the way to 


a consistent treatment of the low energy regime of the 
theory within perturbation theory. 

The general idea was tested in [ 6 l|, where the 
quenched gluon and ghost propagators were computed 
at one-loop order. The comparison with lattice data is 
very satisfactory, with a discrepancy of at most 10 % in 
the whole range of energy. The analysis of the 3-point 
correlation functions (three gluons and ghost-gluon) was 
performed under the same line of investigation in (^ . 
The comparison with lattice data is also very good al¬ 
though the lattice data are more noisy in this case. 

More recently, we have considered the influence of dy¬ 
namical quarks on the theory, which are governed by the 
Lagrangian: 

Nf 

i=l 

where 7 ^ are euclidean Dirac matrices satisfying 

the flavor index i runs over the Nf 
quark flavors^ and the covariant derivative of the quark 
field reads 

- igA’^^F'ip. (3) 

In we have included the one-loop contribution of the 
quarks to the previously computed gluon propagator and 
have computed the quark propagator to the same or¬ 
der. Again, the gluon propagator compares very well 
with the lattice results.^ The quark propagator can be 
decomposed in two independent structures in spinor in¬ 
dices, often called the mass (proportional to the identity 
in spinor indices) and field renormalization (proportional 
to a Dirac matrix). The one-loop results for the quark 
mass is in good agreement with the lattice data. In par¬ 
ticular, it reproduces the strong enhancement at energies 
smaller than roughly 2 GeV, which is reminiscent of the 
chiral symmetry breaking However, we do not expect 
that the present simple one-loop calculation can break 
spontaneously the chiral symmetry. On the other hand, 
the field renormalization of the quark is not reproduced 
correctly, even at a qualitative level. This can be at¬ 
tributed to the fact that the one-loop contribution to 
this quantity is unusually small (^ . Indeed, when the 
gluon mass is set to zero, the field renormalization has no 
contribution at one loop but it does have contributions 
at two-loops, that, in consequence, become dominant. 


^ In this article, we mainly consider degenerate quark masses. In 
this case, the sum over flavors can be replaced by Nf. 

^ Note that there is no direct influence of the quarks to the ghost 
propagator at one-loop. However, the presence of the dynamic 
quarks modify the /9 functions, which gives an indirect effect 
on the ghost propagator. Again, the comparison of the ghost 
propagator with lattice data is very satisfactory. 

® A similar behavior can be seen in the present study in the 
quenched approximation, see Fig |4] below. 
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FIG. 1: One-loop Feynman diagrams for the quark-gluon ver¬ 
tex. 


This model has also been used in another line of re¬ 
search in a recent series of articles: the extension of Yang- 
Mills theory at finite temperature was considered within 
the same framework with the goal of studying the de¬ 
confinement transition It was found that a sim¬ 

ple one-loop calculation can reproduce the known phe¬ 
nomenology [second order transition for the SU(2) group 
and first order transition for SU(3)]. A next-to lead¬ 
ing calculation in the SU(2) case leads to a quantitative 
agreement for the critical temperature. 


III. ONE-LOOP CALCULATION 


The determination of the quark-gluon vertex at one 
loop requires to compute the two diagrams depicted in 
Fig- HI The result can be decomposed in twelve indepen¬ 
dent tensorial structures. We follow here the convention 
of [ 2 ^ and express the quark-gluon vertex as 

=t°'T^{p,r,k) (4) 

with 


^tj.(p,r,k) = -ig 



( 5 ) 


We give in Table U the tensorial structures and 
together with some informations concerning the scalar 
functions Xi and r^. 

The calculations are performed by following the same 
procedure as described in (^.1^. The Feynman rules are 
implemented in a Mathematica procedure and the inte¬ 
grals over momentum written in terms of the Passarino- 
Veltman integrals (^ : 


d‘^q 


(27r)'^ ( 7 ^ -I- 
d‘^q 


A(m) = J 

Bo{mi,m2) = J 

Co{pi,P2,mi,m2,m3) = 

f d‘^q 


{2ttY -I- mf][(p -I- qY + ml 


( 6 ) 

( 7 ) 

( 8 ) 


{2t:Y + "if][(pi + qY + mlWipi +P2 + qY + "^1] 


Our expressions for the different scalar functions ex¬ 
pressed in terms of these integrals are available in the 
supplemental material (b^ . In actual calculations, it is 
convenient to use the Feynman trick to perform the in¬ 
tegrals. In order to reduce the number of Feynman pa¬ 
rameters, we use the identity: 

I _ 

{q^ + m\){q^ + ml) 

ml — ml \q‘^ + m: 

as many times as possible to reduce the number of prop¬ 
agators appearing in the diagrams. The price to pay is 
that there appears differences of masses in the denomi¬ 
nator. For general momentum configurations, after per¬ 
forming the integrals over the internal momentum, we 
are left with a double integral over Feynman parameters. 
However, for integer space dimension and for some mo¬ 
mentum configurations that are considered in next sec¬ 
tion, all integrals can be realized analytically. We thus 
obtain analytic expressions for the one-loop vertex for 
those configurations. 

We have made several checks of our one-loop expres¬ 
sions. We have compared our results with those of (T^ 
which were obtained in the limit of vanishing gluon mass. 
Our check consists in comparing numerically our results 
for vanishing gluon mass with those of [^. This was 
done for several momentum configurations and we find 
a perfect agreement. Therefore, in the regime where 
the external momenta are large compared with the gluon 
mass, our results reproduce those of standard perturba¬ 
tion theory (including anomalous dimensions and all an¬ 
gular dependences). As a consequence, we do not study 
this regime further below. 

Another check comes from the behavior of the coupling 
constants in the chiral limit (limit of vanishing quark 
mass). In this limit, we easily derive the Ward identity 
associated with the chiral symmetry: 

r^75 +75r^ = 0, (10) 

where 75 = 70717273 - By using the anticommutation 
relations of Dirac Matrices, it is easy to check that, in 
the chiral limit, the structures ^ 4 ^, Ti^, 
and are not compatible with this Ward identity. We 
therefore conclude that the associated coupling constants 
must vanish in the chiral limit. We have explicitly verified 
that our one-loop expressions are compatible with this 
result. 

We regularize the theory by using the Infrared-Safe 
(IS) scheme 0- As usual, we introduce renormalized 
quantities, which are related to the bare ones by renor¬ 
malization factors: 


Ab = \fZAA, 

cb = 'J'ZcC, 

Cb = 'J~Zf.c, 


Yb = 

Yb = \^^Y, 


(11) 

to 

II 

2 ry 2 

mB = Zm^m , 

Mb = ZmM. 

(12) 
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Coupl. cte. 

Tensorial structure 

Dim. of cte 

Small M 

Ai 

Llg = 'Tg 

(GeV)'' 


^2 

L 2 g = -{f-f){p-r)g 

(GeV)-= 

oIm°) 

As 

Lig = -i{p - r)g 

(GeV)-i 

oIm^) 

A 4 

L4g = —iogu{p — r)u 

(GeV)-i 

oIm^) 

Tl 

Tig = iikgT^K - rgk'‘‘) 

(GeV)-^ 

oIm^) 

r2 

T2g = {kgTvk,, - rgk'^){f - f) 

(GeV)-'^ 

oIm°) 

T3 

Tig = ^kg - k'^^g 

(GeV) 

o\m°) 

T4 

TA/j, — 'i'\k o'gi_i/{jp 

(GeV)-® 

o\m^) 

n 

Tig = ia gvki, 

(GeV)-i 

o\m^) 

T6 

T&g = ${p - r)g - k^{p - r)^'yg 

(GeV)-= 

o\m°) 

T7 

Trg = -^kx{p - r)x\{f - fhg - {p - r)^] - i{p - r)ga^xr,,kx 

(GeV)-® 

o\m^) 

rs 

Tag = -jgcyi.xr^.kx + - fkg 

(GeV)-^ 

o\m°) 


TABLE I: We give the different tensorial structures on which the vertex was decomposed and the name of the associated 
(scalar) coupling constant. The last line describes the behavior of the coupling constant in the chiral limit. 


The different renormalization factors are defined by im¬ 
posing: 




= fJ-) = M 


^c-co{P = p)= 

r 

ZaZcZ^2 = 1, 

ZgV ZaZc = 1- 


(13) 

(14) 

(15) 

(16) 
(17) 


The last two conditions are consequences of non¬ 
renormalization theorems which relate the diver gen t 
parts of the different renormalization factors (lo - l^ . 
Here, these relations are imposed also for the finite parts. 
The last renormalization condition corresponds to the 
Taylor scheme and is obtained by defining the coupling 
constant through the ghost-gluon vertex function with 
vanishing momentum for the external ghost. The condi¬ 
tion (1161) can not be expressed directly in terms of vertex 
(or correlation) functions. In particular, the mass param¬ 
eter thus defined is not directly related to the value of the 
gluon propagator at vanishing momentum. Expressions 
for the renormalization factors can be found in 0,111. 


In the following, we use the renormalization-group im¬ 
provement for the quark-gluon vertex. It is obtained by 
integrating the RG equation for the quark-gluon vertex 
function 


ldgdg~\~ldf^2djfi2-\- 


Nf^ndM 


(18) 


= 0 , 


where 


Pg{g,m'^,M) 




lA{g,m^,M) 


dg 

% 

drn^ 


go,ml,Mo 


dg. 


go,ml,Mo 


dlogZA 


dg 
d\ogZc 


dg 
dM 


go,ml,Mo 
go,ml,Mo ’ 


dg 


go,ml,Mo 


dlogZ^ 


dg 


go,ml,Mo 


By integrating Eq. (fTSl) we relate the quark-gluon vertex 
at different scales through the equation: 


^tJ.{{Pi},Po,9{po),m'^{po),M{go)) = 

P, gjp), M{g)) (19) 

^ZA{g)z,^{g) 


Here g{g), rn^{g) and Mi[g) are obtained by integrating 
the beta functions with initial conditions given at some 
scale go and: 

logZAig) = ( ^lA {g{g'),m^{g'),M{g')) , 

J go P 

\ogZo{g) = [ ^-fc{g{g'),m^{g'),M{g')) , 

Jgo P 

\ogz^{g)= f ^-/^{g{g'),m^{g'),M{g')) . 

J go P 


Observe that the running quark mass appears explicitely 
in the previous expressions. This leads to a back-reaction 
of this running mass on the vertex functions. As a con¬ 
sequence, the mass takes very different values at /i = 0 
and g = 1 GeV. Even at a perturbative level, this is a 
large effect and, as discussed in the next section, is at the 
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origin of the enhancement of some scalar components of 
the vertex. 

In the IS scheme, for the initial conditions of inter¬ 
est for comparing with the lattice data, the coupling 
constant does not present a Landau pole. In order to 
avoid large logarithms, we evaluate the right-hand-side 
of Eq. (fTOl) with /i of the order of the external momenta 
when they are larger than the mass. In practice, we used 
/i = v^(p2 + k‘^)/2 + m?. It is important to stress 

that contrarily to what happens in the pure glue sector, 
there is an significative dependence on the scheme in the 
present calculation. In fact, it is important to chose the 
renormalization point at a scale that remain fixed when 
all momenta go to zero as done in 0. 

In the following, we compare our results with lattice 
data which were obtained in the quenched approxima¬ 
tion. Consequently, the beta functions for g and m as 
well as the anomalous dimensions for A and c are eval¬ 
uated at Nf = 0. Therefore, the quark mass does not 
back-react on these quantities, but it does back-react on 
the vertex. 


This configuration includes as particular cases the pre¬ 
vious ones but for general an gles only a single tensorial 
structure has been simulated (also in the quenched 
approximation). For the comparison with the simula¬ 
tions, we focus on two values of modulus of the gluon 
momenta and also for the completely symmetric config¬ 
uration where the incoming momentum of the gluon is 
equal in modulus to the quark and antiquark momenta. 
It must be pointed out that the employed data are more 
that ten years old. More recent data with improved lat¬ 
tice parameter could allow to perform a more precise test 
of our curves but they do not exist for the moment. 

In the end of this section we finally discuss the ef¬ 
fect of including dynamical quarks for the present ver¬ 
tex. Even if there are no simulations of this vertex in the 
unquenched case, there are some simulations of 2-point 
functions (t^ that allows to fix the parameters and, in 
consequence, we are able to make a prediction of the 
behaviour of the quark-gluon vertex for the conditions 
that are present in those simulations. These could be 
compared to future unquenched simulations of the quark- 
gluon vertex. 


IV. RESULTS 


The results presented in the previous section were ob¬ 
tained for an arbitrary momentum configuration and for 
all tensor structures. There have been lattice studies of 
the quark gluon vertex for general kinematics for 
a particular tensorial structure in the quenched approxi¬ 
mation. Our results are in general qualitative agreement 
with those of (2^ at the level of precision of the data. 

For some particular momentum configurations, lattice 
data were obtained with larger statistics and therefore 
smaller error bars. In the following, we concentrate on 
these configurations which give more stringent tests of 
the quality of our results. We perform a comparison with 
lattice data in three momentum configurations. The first 
one corresponds to configurations where the gluon has 
zero momentum with arbitrary momentum of the quark. 
In this configuration many tensorial structures were ob¬ 
tained in the quenched approximation and we com¬ 
pare below our results for all available functions. The sec¬ 
ond configuration corresponds to equal momenta for the 
qnark and antiqnark (and, accordingly, glnon momentum 
equal to minus the double of them). Again, many tenso¬ 
rial strnctures were obtained in the quenched approxima¬ 
tion for this configuration (2^ and we study all of them 
below. On top of these two families of configurations, 
we consider a more general configuration with equal mo¬ 
mentum in modulus of quarks and anti-quarks but for 
various possible angles between them (or, equivalently, 
different valnes of the modulus of the gluon momentum). 


A. Vanishing gluon momentum 


Let us first consider the case of zero gluon momen¬ 
tum. In this case, the tensor structure simplifies to the 
following form 


r/.(p,-p,0) = -ig [Ai(p2)7^ - ^\ 2 (p^)fp^ - 2a3p^] . 

( 20 ) 

Moreover, all integrals appearing in the vertex functions 
can be performed analytically. By introducing 


m = 




(Yi — ^ — {M^ — 

{Yi + p^) 2 — (M^ — m^)^ 


(jn? — 

{m? — -I- p2)^ — 


and 

Yi = \J (m2 - M2 -h p2)^ -h 4MV, 
we obtain in d = 4 — e 
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Ai(p^) = 1 + 


327r2 
.2 


- log(m) 


+ 192^2SiVp4 log (^) {-6^® - 1) + [M^ (lOiV^ - 9) - (MiV^ + 3) /] 


-3m'^ - 14M2p2 + 3 + 2) + 3m2 (M® - 7MV + 15MV - /) 

-Ar2(M2+p2)4| ff 


+ 


{mV>"i [I2m^ (iV^ - 1) + (6 - STV^) + mV (iSiV^ + 6) 


3847r2m4iVp4Yi 

-6L2 {2rn^ - 5 M 2 + p^) - MV p^) - ^2^ 


+ Li 


m^ + 2m^ (p2 - M^) + (M^ + p^)^ [ 6 m® - l) + (3 - 4iV2) 


+ m^ (STV^ - 3)p^) -MV^ (V - 3 ) + 6m VV + + 3) p'‘ - (mVp^)^ } , (21) 


and 


A2(p") = 


9 


3847r2m^iVp4 


5 


12 m^ _ 1) _ ^2 (- 4^2 _ 3^ (- 2^42 _ ^ 2 ^ _ 2^2 ^^2 p 2 y 

■ log { 12 m® (iV^ - 1 ) - m® (lOiV^ _ 9 ) (^2M^ - p^) + 6 m'^ (MV^ + p^) 


3847r2m4iVp® 

+ m^ [-6M® - 9MV + 6M2 (SiV^ - l) p'^ + 9p®] 

+2N^ (M® + 4M®p2 - 9MV + 19MV + /)} 

- ^ 0 ^ ~ 

+ 768,^m4JVp8 ^ 1 P" - ”’1 

+ m®p2 [2^//"^ (131V2 - 9) + SM^ (l - 27V2^ p^ + (lOiV^ - 9) p'‘] 

- m^ [150M®1V2 0^6 ^^2 + 1) p2 ^4 ^^q ^2 _ 3) p4 2M^ {2N^ - 3) p® + 3p®] 

+ m^ [300M1°7V2 + 2M® (3 - 1491V2) p^ + M® (47V2 + 15) p'^ + 9MV - (41V2 3^ ^8 

- {2N^ + 3) pio] + 150M®7V2V - 2A^^ {M^ + P^f (75M® + + ZM'^p'^ + 3MV +/)} 

L2 (m^ + 5M^N^) (m^ - + p^) , 


1287r2m4iVp2 


= 64,CV “■® (^) ' 

-p2(m2 - M^) (m^ - 3M^N^) - TVV] - (TV^ - l) p^ 

- ^m^ (iV^ - 1) p2 (m^ + M^) + {rr? - M'^f 

- p^Lz (m^ - mV p2) [m^ + iV^ (p2 _ 3^2^j | ^ 

I- 


4 2 

’ m 


When performing lattice simulations, these functions 
are obtained from the simulated vertex r^(p, —p, 0 ) by 
projecting in the various tensor structures in the follow¬ 


ing way 


[illiii: 


Ai(p^) = ^Tr7^r^(p, -p,0) 


( 22 ) 


(23) 


P(i=0, Pi-#0 for 
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•^2(p^) = -p,o)] + Aib^)^ 

M 

(24) 


^3(P^) = ^Xl^^A*^Re[Trr^(p,-p,0)] 


where, in these expressions, no implicit sum over re¬ 
peated indices is meant. It is interesting to note that, 
in order to extract the function A 2 (p^), a contribution of 
Ai(p^) must be added. This issue will be discussed in 
more detail below. 

To compare our results with the lattice simulation of 
we renormalize the model in the IS scheme and im¬ 
plement the corresponding renormalization group as dis¬ 
cussed in Sect. M- It is important to stress that the sim¬ 
ulations of have been done in the quenched approx¬ 
imation. Accordingly, we use the beta-functions of the 
pure Yang-Mills theory derived in for the gluon 

mass and coupling constant. This implies that the cou¬ 
pling constant and the gluon mass can be fixed once the 
gluon and ghost 2-point functions have been fitted. The 
only remaining free parameter is the valence quark mass 
(and the overall normalization factor). To determine the 
best fit parameters for the ghost and gluon propagators, 
we consider the errors functions: 


Xaa — 


1 

4N 

1 

m 






Xcc — (Mo) + Jit. (Pi)) {Jlt.{pi) Jth.{pi)) 

(25) 


We represent these functions in Fig. [5J The best fits are 
obtained for g = 4.2 and mass m = 0.44 GeV at the scale 
1 GeV which compare well with those found in [^. As 
in that reference, this value has been compared to other 
values of the coupling by taking into account the RG 
running, finding a good agreement with best estimates 
up to a 10 % error. 

The valence quark mass has been fixed in order to fit 
properly the quenched quark propagator. More precisely 
we choose the quark mass in order for the constituent 
quark mass to match the lattice value. Doing so, one 
obtains a mass function (the scalar part of the quark 
propagator, se e [6411 1 that reasonably agree with the lat¬ 
tice data from as can be seen in FigU) One see that 
this curve mimick the enhancement of the mass in the 
infrared but not the spontaneous chiral symmetry break¬ 
ing. Indeed, the constituent quark mass vanishes when 
the current mass at 1 GeV tends to zero. The corre¬ 
sponding value of the quark mass is M = 0.2 GeV at the 
scale p = IGeV, and it gives a constituent quark mass 
M{p = 0) = 0.42 GeV. These values tend to overesti¬ 
mate the mass function M (p) in the UV by approxima- 
tively 30% (for example, at 3 GeV we obtain 0.13 GeV 


6.0 

5.5 
5.0 

4.5 
4.0 

3.5 
3.0 

0.30 0.35 0.40 0.45 0.50 0.55 0.60 


5.0 

4.5 

60 

4.0 

3.5 

0.30 0.35 0.40 0.45 0.50 0.55 0.60 

m (GeV) 

FIG. 2: Contour levels for the error functions xaa for the 
gluon propagator (upper pannel, for 19%, 22% and 24%) and 
Xcc for the ghost propagator (lower pannel, for 4%, 7% and 
10 %). 


m (GeV) 




instead of 0.10 GeV) and to underestimate the propa¬ 
gator in the IR (below 1 GeV) by a maximum value of 
approximatively 18%. Having said that, the overall be¬ 
havior is approximatively correct. It must be pointed 
out when fixing the parameters in order to compare both 
to the quark mass function and to the vertex functions, 
a practical difficulty takes place. The simulations have 
been done for both quantities with different mass param¬ 
eters. For the quark mass function many bare masses 
were employed. In the vertex case, the employed values 
for the bare valence quark mass were Mq = 115 MeV 
and Mq = 60 MeV [^. However the authors claim 
that the vertex functions are almost insensitive to the 
bare quark mass and employ the data corresponding to 
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Mo = 115 MeV. In order to have a single set of data we 
choose the bare valence quark nearest to this value that 
has been used for propagators which is Mg = 75 MeV. 
This introduces a systematic error but that do not seem 
to be large because for relatively small masses like those 

the mass function do not seem to depend on Mg too 
much (see [23 and Fig.lH). 

One can also analyze the dependence on the valence 
quark mass of the vertex by studying the error level for 
the function Ai(p^) defined as 

V = -fr.)- 

We also consider the analogous error level for the func¬ 
tion A 3 (p ^). For this last function the total error may be 
misleading as discussed below because it average some 
momenta where the the curve is essentially within er¬ 
ror bars with the lowest momenta where errors are much 
larger. ® They are both shown in Fig. [31 We ob¬ 
serve that, even if the optimum value of the masses do 
not coincide when fitting Ai(p^) and Xsip^), there is a 
region of parameters where both errors are reasonably 
small (that is less than 10%). Moreover, the value of the 
gluon mass extracted from the ghost and gluon propaga¬ 
tors (m = 0.44 GeV) is within an error region with less 
than 20%. It is important to stress that once the param¬ 
eters g, m and M have been fixed, they are used for the 
calculation of the quark-gluon vertex for all momentum 
configurations. Accordingly, they become pure predic¬ 
tions without free parameters of the present calculation. 

The resulting curves for the functions Ai(p^), X 2 {p^) 
and Xz{p^) can bee seen in Fig. [5] We find a good 
agreement for the function Ai(p^). All the point are be¬ 
low a 16% of error except the first point {p = 140 MeV) 
where a larger error is observed (27%). Let us point out, 
however that our curve is systematically below the lattice 
one. The behavior for XgijP') is qualitatively reproduced, 
with a strong increase (in absolute value) at momenta of 
the order of 1 GeV This last curve depends more strongly 
on the quark mass. This can be understood by the ob¬ 
servation that A 3 tends to zero in the chiral limit, see 
Table m The function is very well reproduced (essentially 
within error bars) except for the two lowest momenta 
{p = 140 and 400 MeV) where the errors are larger (37% 
and 33% respectively when comparing to the lattice cen¬ 
tral value). Let us, however, point out that the statistical 
error bars of these lattice points are large ( 22 % in both 
points). 


^ These values for the mass are not small with respect to actual 
physical values of the masses for light quarks. However, they 
are relatively small with respect to typical QCD scales and, in 
consequence, some quantities are not very sensitive to them. 

® As discussed in detail below, the comparison of our results with 
lattice simulations for the function \ 2 {p^) is problematic. Con¬ 
sequently, we do not use this function for fitting the parameters. 


> 

(U 

O 




0.3 0.4 0.5 0.6 0.7 


m (GeV) 



0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 


m (GeV) 

FIG. 3: Contour levels for the quantity xxi (upper panel) and 
XA 3 (lower panel) in the plane of gluon mass vs. quark mass. 
The contour lines correspond to 4%, 4.5% and 5% for x\i a^nd 
10%, 20% and 30% for XA 3 . 



We therefore find that for a unique choice of parame¬ 
ters ( 7 , m and M, we reproduce very well the ghost and 
gluon propagators and reasonably well the mass function 
M(p) and the vertex functions Ai and A 3 , which is a 
highly nontrivial result. 

The situation looks qualitatively different for the func¬ 
tion A 2 (p^), which tends to a constant at zero momentum 
in our approach while it diverges in the lattice data. This 
disagreement was already observed in previous analytical 
treatments of this function [ 2 l|, and was interpreted 
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FIG. 4: Quark mass as a function of the momentum. The 
renormalization group low is initialized at 1 GeV with q = 4.2, 
m = 0.44 GeV, and M = 0.2 GeV. Lattice data IT^I are 
marqued with dots and parametrized by their bare masses: 
Mo = 19 MeV (squares), Mo = 37 MeV (triangles), Mo = 
75 MeV (circles). Mo = 187 MeV (diamond). 

as a consequence of the non-perturbative features of the 
theory. We want to propose a completely different expla¬ 
nation for this phenomenon. Let us go back to Eq. (IMl) 
which is used to extract A 2 from the lattice data. We 
observe that we must subtract from the function Ai(p^) 
a quantity 

^2 = - ^ Im Tr [ 7 ^r^ {p, -p, 0)]. (26) 

logs ^ 

which is directly extracted from lattice data. The differ¬ 
ence is then divided by p^. In our calculation, the fact 
that A 2 tends to a constant originates from a compensa¬ 
tion of these two terms in the limit p —>■ 0, the difference 
being of order p^. We can interpret this constant as a 
compensation of large numbers, a situation which is dif¬ 
ficult to treat numerically. Indeed, a small error in one 
of these two terms would lead to a divergence of X 2 {p^) 
when p —^ 0 but this would be certainly an artifact. 

Note that such a divergence in A 2 (p^) would imply a 
nonanalytic behavior of the vertex function at small mo¬ 
menta [see Eq. ([20])] since the p —>■ 0 limit would depend 
on the direction (in momentum space) which is used to 
perform this limit. In our calculations, which are done 
in the framework of the action m-m, the only non¬ 
analytic behaviors observed so far are a consequence of 
ghost loops (see, for example, (6ll - [^ ). Such contribu¬ 
tions are not present at one-loop order in the present 
vertex and most probably their contributions of higher 
order should be extremely small. 

If the disagreement takes it origin in partial compen¬ 
sations in the extraction of the A 2 (p^) function from the 
formula (0^ . the comparison should be better for the 
quantity that is directly extracted from lattice data (that 
is A 2 ) than that for A 2 itself. In Fig. |6| we compare our 
prediction for A 2 to the corresponding lattice values. We 
observe a very good agreement in both expressions: our 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 


p (GeV) 



p (GeV) 

FIG. 5: Scalar functions of the quark-gluon vertex: Ai (top), 
A 2 (middle) and A 3 (bottom) for a vanishing gluon momen¬ 
tum, as a function of the quark momentum. The line are the 
one-loop results with M = 0.2 GeV. The dots correspond to 
the lattice data of [ 2 ^. 


curve is almost within the lattice error bars except for 
the lowest momentum where the error is around the 17% 
with respect to the central value of simulations. We must, 
however point out that the curve is systematically above 
the central values of lattice data. The agreement between 
the A 2 with lattice data suggests that the true origin of 
the disagreement between analytical and lattice data for 
A 2 lies in the difficulty of obtaining a function from differ- 
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p (GeV) 


FIG. 6: Comparison of the quantities Ai and A 2 as a function 
of the momentum [see Eq. [26]. Our one-loop results for Ai 
(full line) are almost degenerate with those for A 2 (dashed 
line). The higher (blue) points correspond to Ai and the lower 
(grey) ones to A 2 , extracted from f^ . In the latter case, our 
estimation of the error bar is probably pessimistic. 

ences of large numbers. Let us point out that the origin 
of the error is certainly not statistical. If it were so, the 
signal obtained from the lattice data for A 2 (p^) would be 
extremely noisy. Most probably the origin of the error is 
more systematical. Now it is important to recall that the 
lattice extraction of these scalar function requires a care¬ 
ful treatment of various lattice artifacts (see 
A mismatch in the treatment of lattice artifact of the 
longitudinal and transverse parts of the vertex could eas¬ 
ily introduce a systematic error that, when divided by 
p^, gives a large spurious contribution in the estimate of 
M{p^)- In any case. Fig. [HI shows clearly that not only 
the component transverse to the external momentum but 
also the longitudinal one are correctly reproduced by the 
present one-loop calculation. 

B. Quark and anti-quark with equal momenta 

In this section we present the results for the kinematic 
configuration given by equal momenta of quark and anti¬ 
quark (r = p) and, accordingly, k = —2p. In this kine¬ 
matics the quark-gluon vertex simplifies considerably 

-2p) = -ig [Ai(p^) 7^ -b 4T3(p2)(^p^ - p^y^) 

-2lT5(p^)cr^^Pi,] 

and the transverse projected vertex (that is 
P^(p)ri,(p,p, —2p))< has the form 

KiP^) - 2iT5{p^)af,^p, 

where we have introduced A^ = Ai — k^T^. 

In spite of this simplification, for this configuration 
of momenta, integrals can not be done analytically and 


some remaining Feynman parameter integrals must be 
performed numerically. Fortunately, having fixed the pa¬ 
rameters from propagators, the calculation must be done 
only once for a given choice of parameters. 

In this kinematic configuration, the scalar functions A^ 
and Ts are extracted from the vertex as 

(Tr7^r^(p,p,-2p)) 

and 

^5 = i H (TrcT^,r^(p,p, -2p)) . 

In this case, no subtraction is needed to extract A 3 and 
Ts from the lattice data. 

In Fig. 13 we compare our one-loop results with the lat¬ 
tice data using the same initial condition of the coupling 
constant and masses as in the previous case, g = 4.2, 
m = 0.44 GeV and for the value of M = 0.2 GeV at 
/io = 1 GeV. We observe that the function A'^ is well re¬ 
produced by our results (mostly within the lattice error 
bars, except for the lowest momentum point where the 
error is of the order of 25%). The other function T 5 is 
seen to be very sensitive to the choice of the quark mass 
at 1 GeV. Again, this is a consequence of the fact that, in 
the chiral limit, this function tends to zero, see Table HI 
We observe that the qualitative behavior for this function 
is reproduced but the agreement is not quantitative. In 
particular, our one-loop results do not reach such large 
values in the deep infrared where for the lowest value 
of momenta the error reach 45%. Moreover, the curve is 
systematically above the lattice curve. This is similar to 
what has already been observed in the DS results of . 


C. Symmetric configuration 

In this subsection we consider the case where the quark 
and antiquark have momenta with the same modulus but 
can have arbitrary angle between them (or, equivalently, 
that one can vary the value of the gluon momentum). It 
is clear that the two previous kinematic configurations 
correspond to particular cases of the more general situ¬ 
ation considered here. However, more general configu¬ 
rations of momenta were studied in the lattice [ 2 ^, but 
only for a particular tensorial structure. More precisely, 
in terms of the tensorial decomposition ([S]), the function 
that is extracted in the lattice is A'^ = Ai — as in 
the previous case. We limit to three particular kinematic 
configurations. In the first two, the modulus of the gluon 
momentum taking values |fc| = 0.277 GeV and |fc| = 0.838 
GeV. The third corresponds to the completely symmetric 
case where the modulus of the three momenta are equal 
= p^ = k^. As before, initial conditions of the renor¬ 
malization group equations correspond again to g = 4.2, 
m = 0.44 GeV and M = 0.2 GeV &t pq = \ GeV. 


-2t) = -ig 
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p (GeV) 



p (GeV) 

FIG. 7: Scalar function of the quark-gluon vertex: Ai (top) 
and T 5 (bottom) for equal momentum (p) for the quark and 
antiquark. The full line are the one-loop results and the dots 
correspond to the lattice data of [^ . 


The results are presented in Figs. |S] and 13 For 
|fc| = 0.277 GeV the maximal error is of 12%. For 
|fc| = 0.838 GeV, our results are within lattice error- 
bars with the exception of the point of lowest momen¬ 
tum (p = 0.4 GeV) where the error is of 18%. In the 
completely symmetric configuration the curve is almost 
within error bars except for the lowest momentum where 
the error is 25%. At the end of this section we discuss on 
possible origins of these discrepancies. 


D. Unquenched vertex 

We conclude our study by introducing the effect of dy¬ 
namic quarks in our calculations. This influences the 
renormalization factor for the gluon Za and for the gluon 
mass , which both have a contribution from a fermion 
closed loop. On the contrary, the first contribution of dy¬ 
namic quarks to Z,. and to the quark-gluon vertex comes 
at two loops. We conclude that the effect of dynamic 
quarks only influence indirectly our calculation, through 
the modification of the /3 functions for the coupling con¬ 
stant and the gluon mass and for the field renormalization 




FIG. 8: Scalar function of the quark-gluon vertex Ai = Ai — 
for and k = 0.277 GeV (top) and k = 0.838 GeV 

(bottom). The full line are the one loop results and the dots 
the lattice data of ItOI . 



p (GeV) 


FIG. 9: Scalar function of the quark-gluon vertex Ai = Ai — 
k^Ts for the completely symmetric configuration = 

k^. The full line are the one loop results and the dots the 
lattice data of a interpolated to reconstruct the symmetric 
configuration. 
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The unquenched propagators for the quarks, gluons 
and ghosts where studied in , in the framework of the 
Curci-Ferrari action. We use here the parameters that 
were obtained in that article, as best fits when comparing 
with the lattice data in the only case where the mass 
quark can be fixed from the a quark mass function, that is 
Nf = 2+1. The best fit parameters in this case is g = 4.8, 
m = 0.42 GeV and M = 0.08 GeV at po = 1 GeV. The 
corresponding values of the couplings were compared in 
with standard estimates of the coupling giving an 
agreement with an error of the order of 17 %, which is 
similar to the typical overall error in the infrared of the 
one-loop approximation. 

Using these values of the parameter we can now pre¬ 
dict the vertex in the unquenched case. It is interesting 
to note that in the present procedure this is an extremely 
simple calculation once the beta functions have been cal¬ 
culated and the parameters fixed by comparing to the 
2-point functions. This is in contrast to lattice simula¬ 
tions that are extremely costly in the unquenched case. 
This is probably the reason why there are, to our knowl¬ 
edge, no lattice data for this vertex. 

We focus here on the zero gluon momentum but other 
tensorial structures and kinematic configuration could be 
presented as well. We present in Fig. ITUlthe three scalar 
functions Ai, A 2 and A 3 for different number of flavors. 
The presence of dynamic quarks clearly tends to increase 
these functions at low energy. 


E. Possible origin of discrepancies 

The overall agreement of the present one-loop calcu¬ 
lation is reasonably good. Most of the lattice data are 
well reproduced with a good accuracy. It is important to 
stress that most of the calculations are pure predictions 
without free parameters (that are already adjusted via 
2 -point functions) and that they are all obtained by one 
loop calculations. 

However, some discrepancies are present and we ana¬ 
lyze them in the present subsection. Putting aside the 
A 2 (p), the largest discrepancy concerns the tq{p) func¬ 
tion in the equal momenta configuration. Even if in this 
case the order of magnitude of the function and gen¬ 
eral behavior is well reproduced, there is not quantita¬ 
tive agreement. Finally, one observes some quantitative 
disagreement for the function A^ (p) function in the sym¬ 
metric configuration at small momenta. We can say that 
there is a good general agreement with lattice data but 
the quantitative agreement is not as good as for gluon 
and ghost vertices. A similar behavior is observed in the 
study of propagators . 

We give now two possible (probably complementary) 
explanations for this. First, one of the reasons of the suc¬ 
cess of a perturbative analysis in the gluon and ghost sec¬ 
tor is that the relevant couplings (the value of the 3-gluon 
and ghost-gluon couplings) are moderate except, for the 




p (GeV) 



p (GeV) 

FIG. 10: Scalar functions of the quark-gluon vertex in the 
unquenched case. Ai (top), A 2 (middle) and A 3 (bottom) 
for a vanishing gluon momentum, as a function of the quark 
momentum. The curves corresponds to For Nf = 0 (full line) 
and A/=2+l (small dashes). 


3-gluon vertex for very small values of the coupling where 
the presence of the gluon mass, has already suppressed 
fluctuations. However, in the quark sector, even if the 
coupling does not presents a Landau pole, it reaches sig¬ 
nificantly larger values than in the pure-glue sector. One 
measure of this is the value of the coupling Ai (p) or A)^ (p). 
One can define a natural running coupling function for 
this functions as aq^^^-gluon ^ (y^(-p)) 2 ^ghost-gluon('p)^ 
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It is important to stress that even if for p ^ m, all the 
definitions of the running coupling give a universal beta 
function, this is not the case for momenta p ^ m. In 
consequence, the moderate values of the coupling iir a 
certain scheme does not exclude that it may be large 
in some other scheme. In particular, one observes that 
for small values of p the quark-gluon coupling is larger 
by a factor of 2.4 with respect to the ghost-gluon cou¬ 
pling (giving a larger expansion parameter by a factor 
(2.4)^ ~ 5.8. Iir a sense, this is to be expected because it 
is well established (see, for example, i) that the quark- 
gluon coupling that is required for generating Dynami¬ 
cal Chiral Symmetry Breaking (DSCB) is larger that the 
one that is extracted from the ghost-gluon vertex. Now, 
the expansion parameter of the perturbative expansion 
in the present model in the ghost-gluon sector (precisely 
discussed in [^) may reach values of the order of 0.2. 
In consequence, the expansion parameter in the quark- 
gluon sector seem to put the perturbative expansion in 
trouble at very low momenta. This may indicate that 
even in the present model where the ghost-gluon sector 
seem to be treatable perturbatively, some sort of partial 
resummation is needed in the quark sector (or, at least, 
in part of it). These facts also reflects on a certain 
tension when the parameters have to be chosen in order 
to fit all the data at the same time. Even if the overall 
agreement is good with a single set of parameters, a pre¬ 
cise renormalization scheme had to be chosen because at 
very low momenta an important sensitivity is observed. 

This last point is related to a second possible origin of 
discrepancies. In the quark-gluon vertex, the quantities 
that are not sensitive to the DCSB are much better es¬ 
timated that those that are sensitive to it. The function 
Ai(p), gives better agreement than the Xz{p) or the T^{p) 
(which are sensitive to the DSBC). This manifests also 
in the fact that those quantities seem to have a larger 
sensitivity on the renormalization conditions. Now, it is 
not clear in what measure we have properly taken into 
account the role of the DSBC. Even if RG effect gives an 
important enhancement of the quark mass in the infrared 
that quantitatively reasonably reproduces the mass func¬ 
tion measured in the lattice, it is unclear at the moment 
if a proper treatment DSBC will require to use more so¬ 
phisticated techniques as the introduction of composite 
fields for the chiral condensate. This may be at the origin 
of the lower quality of quantities that are sensitive to the 
DSBC. 


V. CONCLUSION 

In this article, we have presented the calculation of the 
quark-gluon vertex at one loop in the massive (Curci- 
Ferrari) extension of the Landau gauge. The results 
were decomposed on the twelve independent tensorial 
structures. We implement the renormalization-group im¬ 
provement so as to avoid large logarithms. We find that 
these analytic results compare qualitatively well with the 
quenched lattice data. In one case (A 2 for a vanishing 
gluon mass), the difference between lattice data (that 
predicts a divergence at low momentum) and our (which 
are finite at small momenta) may be attributed to a dif¬ 
ference of large numbers when extracting this function 
from lattice data. 

At the quantitative level, we observe that the quan¬ 
tities that are not forbidden by chiral symmetry are re¬ 
produced at the level of 10% except for momenta of 
the order of 100-200 MeV where larger errors are found, 
as detailed in the results section. The others are much 
more difficult to reproduce and show strong dependence 
on the renormalization-group scheme. This is probably 
to be attributed to the rather crude treatment of chiral 
symmetry breaking in our calculation. Moreover, in or¬ 
der to fit simultaneously both the mass function coming 
from the quark propagator and the vertex function, a cer¬ 
tain tension in the choose of parameters is observed and 
in order for this to be fulfilled it is necessary to choose a 
precise renormalization scheme. It is interesting to point 
out that even if, at the end, the agreement of the present 
one-loop calculation is reasonably good, the precision is 
not as good as in the ghost-gluon sector. This may be 
related to a larger value of the associated coupling and 
some sort of resummation in this sector could be useful. 
This issues will be addressed in the future. 
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